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Potentials between static quarks and antiquarks from a few lowest representations were evaluated 
in numerical simulations of 4-dimensional pure G2 lattice gauge theory at various couplings. The 
obtained potentials are linearly rising at intermediate distances and their string tensions exhibit 
(approximate) Casimir scaling. This result is in accordance with a model of the vacuum of non- 
Abelian gauge theories with a domain structure, in which the (color) magnetic flux randomly 
fluctuates within a domain, but the total flux in each domain is quantized in units of the gauge 
group center. 
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1. Introduction 



The behavior of the static quark-antiquark potential in S\J{N) gauge theory in three ranges of 
interquark distances belongs to basic features of confinement in the case of neglecting dynamical 
quarks. For small separations, the interaction is dominated by gluon exchange and the potential 
is Coulomb-like. Going to larger distances, there is an intermediate region with a linearly-rising 
potential. Numerical simulations in both SU(2) [Q] and SU(3) [§, ^ lattice gauge theories showed 
that string tensions for different representations were approximately proportional to the quadratic 
Casimir in that region. The phenomenon is called Casimir scaling. In the region of color screening, 
i.e. at asymptotic distances, the string tensions depend only on the N-ality of the representation. 

From the point of view of confinement, G2 gauge theory represents an interesting laboratory, 
and attracted due attention in the last years. The reason is triviahty of its center which allows to 
reassess the role of the group center in the confinement mechanism. In G2, for example, a "quark" 
from the fundamental representation can be screened by "gluons" (from the adjoint representation), 
thus the string tension is zero at asymptotic distances. However, the authors of [Q] showed that the 
fundamental potential is still linearly rising at intermediate distances. 

To explain linearity, Casimir scahng and asymptotic A'^-ality dependence in a unified way in 
both S\J{N) and G2 gauge theories, a simple model of the Yang-Mills vacuum was suggested in 
Ref. [01 . The model assumes a domain structure of the vacuum where color magnetic fields fluctu- 
ate randomly and (almost) independently inside every 2D-slice of the vacuum, i.e. each 2D-slice is 
a bunch of small independently fluctuating subregions. On the other hand, the total magnetic flux 
through each domain corresponds to an element of the center. For example, in G2 all domains will 
be of the vacuum type because of center triviality. As a prediction, the model leads to Casimir scal- 
ing at intermediate distances in G2 gauge theory. In this contribution we present the first evidence 
for Casimir scaling, more extensive results have meanwhile appeared in Ref. [^. 



2. Potential for fundamental representation 

Our lattice implementation of G2 gauge theory is based on the formulation outlined in [^]. 
One uses a very efficient parameterization of G2 suggested by The Wilson action, given as 

5 = -f £ ^0 Re Tr [Pioix)] - | £ ^Re Tr [Pij{x)] , (2.1) 

' xj>0 ' xj> j>0 '0 

was used. Here Pvt{x) represents the standard plaquette. We worked on asymmetric lattices of the 
type X (2L). The ratio of the lattice spacings in time (a,) and space (a^) directions was fixed by 
requiring the same total physical length of the lattice in all directions, i.e. ^ = cP^^^ /df^^^ = 2. The 



bare anisotropics t,Q, entering the Wilson action (gj) and leading to = 2, depend on the coupling 



j3 and were determined following the procedure of Ref. [||]. 

The simulations were performed at three values of the coupling j8 € {9.5,9.6,9.7}. For each 
coupling we produced almost 1000 configurations. In order to have overlap of our trial quark- 
antiquark state with the ground state as large as possible, we applied the so-called stout smearing 
procedure of Morningstar and Peardon properly adjusted for G2. We optimized the smearing 
separately for each potential. 
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Figure 1: (Left:) Potentials Vst, Vss, and Vrs. Distances r are measured in units of as- 14^ x 28 lattice, 
/3 = 9.6, = 1.97. (Right:) Potentials V^, for representations {D} as functions of the dimensionless r. 
143 X 28 lattice and j3 = 9.6 

The static potential is expressed by means of the expectation values W{r,t) of Wilson loops as 



V{r,a) 



-lim-lnW(r,f) 

t^oo t 



(2.2) 



where r = r • a is the spatial extent and t = t ■ a is the temporal separation of the Wilson loop. 
Because of various orientations of the Wilson loop on an asymmetric lattice, one can define three 
different potentials - V^t, V^s and Vts- They are not independent, and each potential V{r,a) consists 
of two contributions: the true r-dependent static potential due to the interaction, and the self energy 
contribution (independent of r). The interaction parts satisfy two simple relations (see |8|]), 
allowing us to determine the bare-anisotropy parameter of the action. 

In accordance with the relation ( |2.2D , the potentials (in lattice units) were determined by fitting 
logarithms of the measured loops by linear functions in t: 



lnW;,{r,t)=C^+V;,{r)-t, A G {ss,st,ts}, 



(2.3) 



in an interval of t values, from fmin to t^nax, where \nWx{r,t) is approximately linear. Using the 
smearing procedure, logarithms of Wilson loops show a linear behavior down to f = 1, so the lower 
Umit was chosen tmm = 1. The upper limit of our fits tmax was not fixed so firmly, and was varied 
depending on the ground-state overlap. By means of the standard least-square method, we esti- 
mated statistical errors of the linear fit ( ^j| ) in a fixed (fmim^max) -interval. Systematic uncertainties 
and errors of the procedure were discussed in more detail in Ref. [^]. 

In the second step, the calculated potentials were fitted by a 3-parameter function, including 
the Coulomb and linear terms: 

"A 



+ a^r, X €^ {ss,st,ts}. 



(2.4) 



For illustration, all three potentials are displayed in Fig. [I] (left panel) for j3 = 9.6. We can identify 
easily the range of the linear behavior, and we can also conclude that the potential Vst is most 
suitable for reliable determination of the string tension. This conclusion is supported by results of 
the potential-fitting procedure, summarized in Table |l|. 

Finally, from the numerical values of the string tensions we computed the renormalized aniso- 
tropics ^ = 2.00(7), 1.97(12), and 1.96(16) for j3 = 9.5,9.6, and 9.7, respectively. These values 
confirm that the renormalization was done accurately. 
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p 


ttst 








CCts 




9.5 
9.6 
9.7 


0.118(5) 
0.124(4) 
0.122(5) 


0.1997(19) 
0.1303(18) 
0.0999(21) 


0.249(31) 
0.264(34) 
0.264(36) 


0.399(13) 
0.257(15) 
0.196(15) 


0.194(5) 
0.231(7) 
0.244(8) 


0.2482(22) 
0.1784(30) 
0.1452(33) 



Table 1: Parameters of the static potential between fundamental charges. (14^ x 28 lattice.) 



3. Potentials for higher representations 

Construction of higher representations: Using tensor decompositions of different products 
of representations, traces of higher-representation matrices can be expressed through traces of the 
fundamental- and adjoint-representation matrices, e.g. in the case of the 27 -dimensional represen- 
tation one obtains the relation 

TrD^^'^^ = - l-TrD^-TrD^ + (TrD^)^ (3.1) 

The adjoint-representation matrix is constructed from the fundamental one: 

Bf,,{g) = 2 Tr [n^igy ta J}^{g) th] , (3.2) 

where ta are the generators of the group, based on representation described in [Q]. 

Determination of potentials: We focus our attention only on potentials which can be cal- 
culated most accurately and which exphibit best the linear behavior. On a thermalized, smeared 
configuration in the fundamental representation, we computed a Wilson loop as a product of link 
matrices in the fundamental and adjoint representations, and knowing these two traces we calcu- 
lated also traces of the Wilson loop in higher representations. Finally, we computed their expecta- 
tion values. 

To obtain values of potentials, we use again the linear fit (pj|). The potentials in higher repre- 
sentations show similar behavior to the case of the fundamental representation - a linear rise in the 
same range of distances, but with different slope (see Fig. [I], right panel). 



Finally, we fitted the computed potentials again by the 3-parameter fit ( |2.4| ). The resulting 
ratios of the string tensions are summarized in Table ^. The calculated ratios are very close to pre- 
dictions based on Casimir scaling, and we observe an improvement of the agreement with Casimir 
scaling by increasing the coupling j8 . Fig. ^ shows ratios of different-representation potentials to 
the fundamental one at fixed r values. For each r, the ratios barely deviate from the Casimir-scaling 
prediction. 





Oa/Of 


5{27}/5f 


^{64} /^F 


5{77}/5f 




9.5 


1.88(4) 


2.15(5) 


3.1(1) 






9.6 


1.94(4) 


2.24(6) 


3.35(8) 


3.8(2) 


4.6(2) 


9.7 


1.96(6) 


2.28(7) 


3.5(1) 


4.0(2) 


4.9(2) 


Prediction 


2.0 


2.333 


3.5 


4.0 


5.0 



Table 2: Ratios of string tensions Gst for different representations. (14^ x 28 lattice.) 
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Figure 2: Ratios Vv,.{£)}/Vw,f for different representations {D} are shown as functions of the dimension- 
less r. Horizontal lines represent predictions of Casimir scaUng. (14^ x 28 lattice.) 

4. Conclusions 

In the G2 lattice gauge theory, we calculated string tensions for static potentials between color 
charges from the six lowest representations of G2 on asymmetric x (2L) lattices. Ratios of the 
string tensions in the interval of intermediate distances exhibit Casimir scaling with very good 
accuracy. The deviations from predictions based on values of quadratic Casimirs were estimated - 
including both statistical and systematic errors - to be at most 10-15%. The results, together with 
the evidence for Casimir scaling in SU(2) and SU(3), support the model of magnetically disordered 
Yang-Mills vacuum with a domain structure. 
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